Quantum logic via the exchange blockade in ultracold collisions 
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A nuclear spin can act as a quantum switch that turns on or ofT ultracold collisions between 
atoms even when there is neither interaction between nuclear spins nor between the nuclear and 
electron spins. This "exchange blockade" is a new mechanism for implementing quantum logic gates 
that arises from the symmetry of composite identical particles, rather than direct coupling between 
qubits. We study the implementation of the entangling VSWAP gate based on this mechanism for 
a model system of two atoms, each with ground electronic configuration ^So, spin 1/2 nuclei, and 
trapped in optical tweezers. We evaluate a proof-of-principle protocol based on adiabatic evolution of 
a one-dimensional double Gaussian well, calculating fidelities of operation as a function of interaction 
strength, gate time, and temperature. 
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Implementations of quantum processors require coher- 
ent quantum logic gates that induce entangling opera- 
tions between qubits [l|. In most cases, the same physical 
effect that gives rise to coherent couplings in the logical 
basis states leads to decoherence between them. For ex- 
ample, in an ion trap, the strong Coulomb interaction 
provides the quantum bus for two-qubit logic 0], but 
charged ions can be heated by fluctuating patch poten- 
tials in trap electrodes [1]. In the neutral atom platform, 
strong spin dependent collisions can lead to an entangling 
gate [J, 0, Q but also to spin relaxation Q . 

In this article, we propose a new scheme for quantum 
logic based on an "exchange blockade" arising solely from 
the symmetry of identical composite particles rather than 
from differential coupling strengths. We consider a hy- 
brid approach based on NMR and ultracold collisions of 
trapped neutral atoms. Elements with two valence elec- 
trons in a closed shell configuration ^Sq, such as those 
in group-II, possess no electron spin and thus no hy- 
perfine interaction with nuclear spin when in the elec- 
tronic ground state. Nuclear spin, nonetheless, strongly 
affects the interaction between these atoms solely due to 
quantum statisitics. In a fermionic species, the dom- 
inant ultracold s-wave collisions are forbidden by the 
Pauli exclusion principle when the nuclei are in a sym- 
metric spin state, and are allowed only for the nuclear- 
spin antisymmetric states. The single-channel cold colli- 
sions, governed by electronic interactions between atoms, 
thus induce an exchange interaction between nuclei, even 
though there is no direct interaction between nuclei nor 
between electrons and nuclei. The ability of nuclear spins 
to act as a switch that turns on and off electron interac- 
tions provides a mechanism by which quantum informa- 
tion can be protected from the environment while simul- 
taneously qubits are strongly coupled to one another in 
a manner depending on the logical basis states. 
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We take group-II-like elements with isotopes of spin 
7=1/2 nuclei to encode qubits, e.g. ^^^Yb, with logical 
basis states defined in the usual way, |0) = | t), |1) = | i)- 
Single qubit rotations can be achieved via rf-pulscs in a 
strong bias magnetic field, as in NMR. Interactions be- 
tween atoms that have nuclei spin polarized in the same 
direction are blockaded (to the degree that p-wave col- 
lisions are negligible). Atoms undergo an exchange in- 
teraction as the nuclear-spin singlet component receives 
an elastic collisional phase shift relative to the nonin- 
teracting triplet component, resulting in the well known 
VSWAP, a universal entangling gate pj]. The exchange 
interaction was previously studied in the seminal work 
of Loss and DiVincenzo in the context of electrons in a 
double quantum dot Q. In our realization, the compos- 
ite nature of the atom allows one to separate the degree 
of freedom which stores the quantum information (the 
nucleus) from that which provides physical coupling (the 
electrons). 

Coherent collisions require well-localized atomic wave 
packets whose motion is highly controlled. We evalu- 
ate the performance of two-qubit entangling gates using 
dipole traps formed by tightly focused optical tweezers, 
recently used to trap individual atoms @. Dorner et 
al. have studied ultra-cold collisions of bosonic alkalis 
in such tweezers [l3]- Similar gates have been studied 
theoretically in a lattice of double- well potentials formed 
by multiple frequencies fX)\ and more recently, experi- 
ments explore the control of atoms in double- well lattices 
formed in 2D through overlapping lattices of orthogonal 
polarizations fl2j . We revisit the problem of implement- 
ing gates in double- well potentials, here in the context of 
the exchange blockade. 

We take a simple one-dimensional model assuming 
tight confinement in directions transverse to the direction 
of the collision. We model the two-tweezer dipole trap as 
a symmetric double well of Gaussian shaped potentials, 
each of depth Vq and rms width cr, separated by distance 
d. Single atom orbitals are the eigenstates of the double 
well, which are then filled with two atoms for a two- 
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qubit system, in analogy to the two-electrons occupying 
the molecular orbitals of a diatomic hydrogen molecule. 
The two atomic qubits interact via a contact term arising 
solely from s-wave collisions, V^nt(a;i — X2) = g5{xi — X2)- 
The coupling constant, 5, follows from a quasi-lD ap- 
proximation to the scattering process, ignoring renormal- 
ization effects 13]. Throughout, we choose harmonic 
oscillator units such that energy is measured in units 
of huj, where uj — ^JVoJJjria^ is the oscillation fre- 
quency along X, and lengths are measured in units of 
Xq — ^Jh/{muo). The dimensionless interaction strength 
is then g = 2{as/xQ){ijj]^luj), where Og is the s-wave 
scattering length and u)±_ is the oscillation frequency 
in the transverse dimension. We take a well depth of 
Vo = lOSo;, or hw — 20Ei^/{ka)^, where Efi is the recoil 
energy. Though in principle deeper wells and larger oscil- 
lation frequencies are possible , high fidelity operation 
requires optimal control of that depth and an asymmetric 
"tilt" of the well [10, 12], not under consideration here. 

To implement a gate, we envision separated noninter- 
acting atoms, each in its own tweezer, that are brought 
together into the same well as d ^ where cold collisions 
occur and then are separated again into individual traps. 
As the wave packets overlap, atoms lose their identity, 
and proper Fermi symmetrization of the state is neces- 
sary. We define our two-qubit logical basis such that the 
first qubit is localized in the left (L) well and the second 
in the right (i?), 

|0,0) - fl^fljvac) = ^-ix,,X2)\ TT), (la) 
|0,1> = /It/LHc) (lb) 

= ^{■^-{xi,X2)\xt) +-^+{xi,X2)\xs)) , 

= ^{<S-{xi,X2)\xt) - ■^+{xi,X2)\xs)) , 

= ^ *-(2^i.^2)l ii). (Id) 

Here we have used both first and second quan- 
tized notation, where P are fermionic creation 
operators for localized atoms, '^±{xi,X2) — 
{'>Jjl{xi)iPr{x2) ±^Rixi)ipLix2)) /V2 are two-atom 

orbitals,^ and ]xt(5)) - (1 T)l i) ± 1 i)l T)) /V2 are the 
triplet (singlet) nuclear spin states with zero projection 
of angular momentum. The exact motional state is 
irrelevant to the encoded quantum information in the 
nuclear spin, but may affect the nature of the collision 
and fidelity of the gate. We begin by taking the atoms 
to be cooled to their motional ground state and later 
treat the effects of finite temperature. 

As seen in the logical basis states, Eq. ((TJ, a SWAP 
operation, ]0, 1) <^ 11,0), follows from a sign change be- 
tween singlet and triplet components, or a tt rotation 
in this subspace; the entangling VSWAP is a tt/2 ro- 
tation. As a first attempt to implement this gate, we 
consider adiabatic evolution, though certainly this is not 



optimal for speed. In order to better understand this, 
consider the case of noninteracting atoms. Single par- 
ticle eigenstates are parity eigenstates, which for large 
separation, become doubly degenerate pairs of states cor- 
responding to gerade and ungerade combinations of the 
left /right localized wave packets in the separate tweezers 
denoted \gn{un)) — i\Ln) ± \Rn)) /V2 for the n*'' vibra- 
tional eigenstates. For zero separation, \gn) \v2n), 
\un) \v2n+i) , whcrc \vn) IS a vibration eigenstate of 
the unified well. At large separation and in the absence of 
interactions, there is a four-fold degenerate ground-state 
manifold of two-particle orbitals. These are described 
by the spatial orbitals with one atom per well as de- 
fined in the logical basis \E'±, and those with two-atoms 
per lattice site (associated with spin singlets), ]^'c) = 
(]L)]L) + \R)\R))V2 and ]«-,) = {\L)\L) - \R)\R))V2. 
The state ]^d) has odd parity under spatial inversion 
and must be antisymmetric under spin exchange, and 
thus does not couple to the even-parity spin singlet in 
the logical basis. 

In the presence of interactions, degeneracy between 
singly and doubly occupied states is broken. Localized 
states become eigenstates for separated traps as in the 
Mott insulator phase [1, [l3]. The degeneracy of states 
]5'+) and ]^'c) for separated wells is broken by the onsite 
interaction. A relative phase then develops between sin- 
glet and triplet c omponents of the logical basis, and the 
desired vSWAP operation is achieved when, 

=vr(l/2 + n), (2) 

where (/)± = dt' E±{t')/h is the adiabatic dynamical 
phase. 

Adiabatic energy curves (relative to the noninteract- 
ing ground state singlet adiabatic curve) for states that 
have no more that one quantum of vibration in sepa- 
rated wells are shown in Fig. 1, for different interac- 
tion strengths. Only even-parity singlets and odd-parity 
triplets are shown since these are the symmetries of the 
logical basis states. The triplet states are unchanged by 
interactions and singlets are found by diagonalizing the s- 
wave contact interaction in the noninteracting basis. The 
basis orbitals are calculated using a DVR method 15f] to 
find single-atom bound states in the double well. These 
are occupied with two atoms in accordance with Fermi 
symmetry. 

The behavior of the adiabatic curves as a function 
of interaction strength elucidates the possible operating 
points of a high-fidelity entangling interaction. For weak 
interactions, g <C 1, the very small energy gap between 
]4'+) and ]4'c) limits the ability to maintain adiabaticity 
as the traps are merged. In this regime, the tunneling 
energy dominates over the interaction energy before the 
wells are fully merged, and diabatic transitions arise as 
the ground state wave function changes from right/left 
localized to delocalized across the double well. In prin- 
ciple, tunneling can be suppressed through an asymmet- 
ric well "tilt" that provides stronger symmetry breaking, 
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FIG. 1: Adiabatic energy curves relative to the noninteract- 
ing ground state as a function of trap separation for the first 
two triplet (blue dotted lines) and the first three singlet states 
(red solid lines) and for different values of interaction coupling 
strength g: (a) g = 0, (b) <? = 0.2, (c) 5 = 1, (d) = fO. For 
large separation, these asymptote to single well energy levels 
with or 1 unit of vibration. All quantities are normalized 
to harmonic oscillator units (see text). For large interaction 
strength, (d), singlets and triplets become degenerate corre- 
sponding to the Tonks-Giradeau regime. 



with the capacity to maintain adiabaticity while simul- 
taneously inducing an exchange interaction (lol . [l^ . For 
intermediate interactions, g 1, a large gap is opened 
between l^'-f) and \^ c) , in which case symmetric double 
wells is the preferred geometry for adiabatic operation 
of the gate. For g ^ 1, whereas one might expect the 
most favorable operation of the exchange blockade, para- 
doxically the opposite is true. Singlet and triplet poten- 
tials again become degenerate as one enters the Tonks- 
Girardaeu regime of the one dimensional potential . 
In a contact interaction with infinitely strong coupling 
in one dimension, the particles are excluded from being 
at the same position, regardless of the symmetry of their 
spin state. As a result, exchange interaction between 
spins is reduced. This is illustrated in Fig. 2 which shows 
the gate time in the case of adiabatic evolution as a func- 
tion of interaction strength for different values of n in Eq. 
(PI). The divergence of r as 5 ^ cx) represents the Tonks 
regime. Since the Tonks gas is an artifact of our one- 
dimensional model, this is not a fundamental limitation. 
In future work we will study tight confinement in 3D. 

Though gates times can be calculated based on Eq. 
(121) , there is no guarantee that these times are consistent 
with adiabatic evolution. To evaluate the performance 
of this protocol, we must numerically solve the full time- 
dependent Schrodinger equation. Expanding in the adi- 
abatic basis {|a(t))} and going to the rotating frame of 
the dynamical phase, the probability amplitudes evolve 
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FIG. 2; Gate time under the assumption of adiabatic evolu- 
tion as a function of g for different values of n in Eq. ((2|; 
n = 1,5,10,25,50,100,200,250,300. Note, for n small, the 
short gates times give very low fidelities due to diabatic ef- 
fects in the true evolution (see Fig. 3). The divergence of 
these curves as g becomes large follows from the effect of the 
Tonks-Girardaeu regime. 

according to, 

^" = -EH|l/3)e-^('^^-^°)c^. (3) 

In our numerical integration, we take a finite basis con- 
structed from the span of two-particle orbitals with a 
maximum of two quanta of vibration in the separated 
wells. The trajectory, d{t), is chosen as a (suboptimal) 
linear ramp so that the vSWAP is achieved according to 
the adiabatic dynamic phase evolution. 

We seek the fidelity of operation representing the over- 
lap between final and target states. Since quantum in- 
formation is stored in the nuclear spin, we trace over the 
vibrational excitations. For identical particles, one must 
take care so as not to artificially introduce entropy as- 
sociated with (anti)symmetrization; subsystems must be 
defined by left /right localized atoms, not by particle label 
The 4x4 reduced density operator for the qubits 
is thus, 

n,n' 

(4) 

where L„ and i?„ are left/right localized wavepackets 
with n units of vibration, and x and y denote or 1 
for the logical basis. In general, this density matrix will 
not have unit trace since there is finite probability of 
finding two atoms in the same well. Double occupancy 
represents an infidelity in the operation. The remain- 
ing infidelity arises from an error in reaching the target 
two-qubit spin state. We define this target as the two- 
atom spin state reached by adiabatic evolution for time 
r, where the initial state is the logical basis state |0, 1). 
The other possible initial logical basis states |0,0) and 
|1, 1) evolve to good approximation by the unit operator 
and do not significantly contribute to infidelities when 
overall performance is good. Our fidelity is the overlap 



4 



between the chosen target and the two-atom spin state as 
calculated by the numerical solution under full evolution 
by the time-dependent Schrodinger equation, 

T ^ ('0targot(T)|/5spi„(T)|'(/'targct(T)). (5) 

When the time r is chosen to satisfy Eq. ([2]), plotted in 
Fig 2, this is the fidelity for achieving the desired gate, 

IV'target) = VSWAP|0,1) - ( 1 0, 1) + 1 1 1 , 0) )/ ^2. 

Figure 3 shows the fidelity of adiabatic operation as a 
function of gate time and interaction strength. Though 
high fidelity can be achieved for a wide range of pa- 
rameters, robust operation will require high interactions 
and long time evolution. As an example, for a coupling 
strength of g = 8 and gate time of r = 36(27r/(jj), the fi- 
delity to implement a vSWAP is 0.982. Though not nec- 
essarily contributing to gate infidelity, vibrational heat- 
ing can occur due to nonadiabatic dynamics. We can 
calculate the mean vibration excitation in one atom af- 
ter the operation according to 

(n) = ^n ('^ac|/L„,.:/fl„,,2,P/l„,:./i;„,,>ac). (6) 

For the example above, the interaction will heat each 
atom to have {n) = 0.016 quanta of vibration, on average. 




FIG. 3: Fidelity of adiabatic operation starting in the logical 
basis state |0, 1), according to Eq. ((5} as a function of g and 
gate time r. When r and g are chosen according to Eq. ([2}, 
this corresponds to the fidelity of the \/SWAP gate. 

The analysis above assumed initial preparation of the 
atoms in the vibrational ground state. We can study 



the tolerance to finite temperature in the initial state by 
averaging over a Boltzmann distribution of vibrational 
excitations. The fidelity as a function of (n(0)) falls off 
approximately linearly, T {{n{Q))) ~ —1.78 (n(0)) -I- J^o- 
The particular value of To depends on the parameters, as 
shown in Fig. 3. A necessary condition for high fidelity is 
that initial vibrational excitation satistify {n{0)) < 0.01 
or kT/hw < 0.1. This extreme cooling requirement limits 
the practically of implementing this protocol in the lab. 
Our focus here is on the proof-of-principle protocol, far 
from optimal, and generally very slow. Note, however, 
because the atom's nuclear spin is decoupled from elec- 
tronic states, in principle damping forces can be intro- 
duced which maintain atoms near the ground vibrational 
state without decohering the qubit. Such refrigeration 
can be introduced by laser cooling that does not flip nu- 
clear spin [20| or through immersion in a superfluid where 
sympathetic cooling can occur [l^ . More optimal proto- 
cols will be studied in future work. 

In conclusion, we have presented a new protocol for en- 
tangling qubits based on an exchange blockade between 
identical particles. Though we have studied this pro- 
tocol in the context of identical spin-1/2 fermions, our 
result generalizes for an arbitrary spin, qubits or qudits, 
Bose or Fermi, and is not restricted to elements with elec- 
tron and nuclear spin decoupled. An equivalent evolution 
will arise if, in the absence of symmetrization for iden- 
tical particles, the scattering phase shift is independent 
of logical basis state. For the case of ^Sq elements with 
an arbitrary nuclear spin, we can define 2s -I- 1 logical 
states {|n, m)} which will undergo a vSWAP operation, 
\n, m) |n, m) -t- i|m, n), as a consequence of a cold col- 
lision and the exchange blockade. As another example, 
two ultracold ^^Rb atoms (bosons) will incur a s-wave 
collisional phase shift essentially independent of the mag- 
netic sublevel. If one encodes a qubit in two sublevels, 
|0) = |F, m),|l) = |i^',m'), the logical basis states are 
equivalent to Eq. ([T]), with bose-symmetry substitution, 
\^+) l^*-)- The s-wave collisions now occur for the 
symmetric combination of internal hyperfine levels, but 
not for the antisymmetric combination. An equivalent 
protocol of bringing identical atoms into the same well 
and then separating them will yield a VSWAP [Toj . 
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